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Ph, Abstract 

C^H. The two loop equal mass sunrise graph is considered in the continuous d-dimensional regularisation 






for arbitrary values of the momentum transfer. After recalling the equivalence of the expansions at 
d = 2 and d = 4, the second order differential equation for the scalar Master Integral is expanded 
in (d — 2) and solved by the variation of the constants method of Euler up to first order in (d — 2) 
included. That requires the knowledge of the two independent solutions of the associated homogeneous 
equation, which are found to be related to the complete elliptic integrals of the first kind of suitable 
arguments. The behaviour and expansions of all the solutions at all the singular points of the equation 
are exhaustively discussed and written down explicitly. 

Key words: Feynman diagrams, Multi-loop calculations 
PACS: 11.15.Bt, 12.20.Ds 



1 e-mail: Stef ano.LaportaObo. infn.it 
2 e-mail: Ettore.Remiddi@bo.infn. it 



1 Introduction. 

The two-loop sunrise selfmass graph with equal masses is not yet known analytically, despite its apparent 
simplicity. In this paper we fill that gap and evaluate it, for arbitrary values of the momentum transfer, in 
terms of a suitable set of new functions, whose analytic properties (in particular behaviours and expansions 
at all their potentially singular points) are worked out explicitly in full details. 

We follow the differential equation approach, already introduced in [1] in the arbitrary mass case for 
studying particular values of the momentum transfer, within the usual regularization scheme in d-continuous 
dimensions. We are interested, as usual, in the d — > 4 limit; but it is known from [2] how to relate 
algebraically the coefficients of the expansions at d — 4 and those of the expansion at d — 2, and as the 
expansion at d = 2 turns out to be simpler, we discuss in fact the d — > 2 limit. 

In the equal mass limit the two loop sunrise has two Master Integrals (Mi's) only, which satisfy a system 
of two first order differential equations in the momentum transfer u (we will use often also the Euclidean 
variable z = —u, which is positive when u is spacelike), equivalent to a single second order equation for 
the simpler of the two Mi's. The equation, which is exact in d, when expanded around d — 2 in powers of 
(d — 2) gives rise to a set of chained inhomogeneous equations all having the same homogeneous part; we 
solve them recursively by the variation of the constants method of Euler, which gives formally the solutions 
of the inhomogeneous equations in terms of the solutions of the common homogeneous equation. Only 
when the homogeneous solutions are known Euler's formula does provide an effective analytic evaluation; 
the bulk of the paper consists, indeed, in working out those homogeneous solutions. 

The paper is organized as follow. In Section [3 the differential equations are written; in Section [3] the 
connection between the expansions at d = 4 and at d = 2 is established; in Section 0] the expansion of 
the second order differential equation at d — 2 is worked out. Section discusses how the solutions of 
the homogeneous equation can be obtained; Section presents the expansions in z of the homogeneous 
solutions at all the singular points of the differential equation; Section [3 discusses another set of functions, 
represented in the form of one-dimensional definite integral over a suitable parameter, which solve the 
homogeneous equation between any two nearby singular points (referred to as "interpolating solutions"). 
Section |H1 then finally gives the proper analytic continuation of the homogeneous solutions for any value 
of z, and Section El discusses the properties of the homogeneous equation and of its solutions for special 
conformal transformations of the argument. 

Once the homogeneous solutions are known, Sectional works out explicitly the MI as given by Euler's 
formula for the solution of the inhomogeneous equation at zeroth order in (d — 2), Section^] the corre- 
sponding zeroth order term of the expansion in (d— 4) and Section IT2l the first order term of the expansion 
in (d-2). 

Section El contains the conclusions; the first of the two Appendices deals with the relation between the 
"interpolating solutions" and the complete elliptic integrals of the first kind, the second Appendix with 
some definite integrals occurring in Sections HUIandlT2l 

2 The differential equations for the Mis. 

The 2-loop sunrise graph of Fig. U with arbitrary masses mi, TO2, TO3, is known to possess 4 Master Integrals 
(Mi's) [3], already extensively studied in the literature [4]. The system of the four linear differential 
equations in the squared external momentum p 2 for the Master Integrals (Mis) was written in [1]. The 
differential equations were used for obtaining analytically specific values or behaviours at zero and infinite 
momentum transfer [1] , at pseudothresholds [5] and threshold [6] , as well as for direct numerical integration 
[7] , but a satisfactory knowledge of the analytic expression of the Mi's for arbitrary momentum transfer and 
masses is still missing. In this paper we tackle that problem in the simpler equal masses case. While the 
algebraic burden in the arbitrary mass case will surely be much heavier, there are however indications [8] 
that the approach can be extended to the arbitrary mass case as well. 

We will work as usual within the d-continuous regularization scheme, defining the loop integration 




Figure 1: The 2-loop sunrise graph. 

measure, in agreement with previous work, as 
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At d = 4, it reduces to the standard measure 
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further, with the definition Eq. (|2.1|l . the tadpole reads (exactly in d) 
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When all the masses are equal, the 2-loop sunrise has just two Mis; putting equal to 1 the value of the 
mass, the two Mi's S(d,p 2 ) and Si(d,p 2 ) can be defined as 

S ( d >P ) = / tu2 , TVL2 , ^r/„ u 1 ^2 , Ti . (2- 4 ) 



(fc? + l)(fc| + l)[(p-fci-*a) 2 + l] ' 



2 /■ D d fciD d fc 



2 



^ (d ' P } _ J (fc 2 + l)2(fcl + l)[(p-fc 1 -fc 2 )2 + l] ' (2 ' 5) 

where p is the external momentum. In the following we will use 

p = z = —u , (2-6) 

with z positive (and u negative) when p is spacelike, and the continuation to timelike values (u positive) 
will be performed with the usual replacement 

z = -(u + ie) . (2.7) 

The system of linear differential equations in z satisfied by the 4 Mi's in the arbitrary mass case was 
written in [1] ; in the equal mass limit it goes into the following 2x2 system of linear differential equations 



in z for the two Mi's defined in Ea.s l|2.4l2.5|) 
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The system can be rewritten as a second order differential equation for S(d, z) only 
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and the second MI Si(d, z) can be expressed in terms of S(d, z) and its first derivative, 
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From now on we can take S(d, z) and its first derivative dS{d, z)/ 'dz as the effective Mis of the problem. 



3 The expansions in d. 

In any physical application, the d — > 4 limit is to be taken; in that limit, S(d, z) can be Laurent-expanded 
in (d — 4) around d = 4 as 

S(d,z)= ]T S<">(4,*)(d-4) n , (3.1) 

n=-2 

where, as it is well known, the Laurent-series at d = 4 starts with a leading double pole in (d — 4). In this 
paper we will show how to evaluate analytically the coefficients S^ n '(4:, z) Up to n = 1, i.e. up to the first 
order term in (d — 4) included. 

Let us recall here that by acting on any Feynman integral in d-continuous dimensions with a suitable 
differential operator one can obtain an expression relating the value of that integral evaluated in (d — 2) 
dimensions to the values, evaluated in d dimensions, of a suitable combination of other integrals related to 
the same Feynman graph [2]. In the case of the Master Integral S{d, z), expressing the result in terms of 
the effective MPs (and of tadpoles) the relation reads 
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By differentiating in z the above relation, and using Ea H2.9|) for eliminating the second derivative of S(d, z), 
one obtains a similar relation for dS(d — 2, z)/dz. The two relations can be used to express S(d, z) in terms 
of Sid — 2, z) and dS(d — 2, z)/dz; replacing further d by (2 + d), one finally obtains 
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Quite in general, assume a relation of the form 

L(2 + d)= R(d) (3.4) 

as given; setting d — 2+rj and Laurent-expanding in 77 around d = 2 both the l.h.s., £(4+7/) = ^ n lS n ' (4)ry n , 
and the r.h.s., R{2 + r/) = J2 n ^ (n) (2)7j n , one obtains 
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which obviously implies, at any order n in the Laurent-expansion in 77 

L^(4)=R^(2). 

We put d = 2 + 77 in Eq. (|3.3() and look for the systematical Laurent-expansion of l.h.s. and r.h.s. The 
l.h.s. is nothing but 
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where the coefficients 5^(4, z) are the same as in Eq. fl3.lj t. Within the r.h.s. of Eq. fl3.3|) we write 

S(d,*) = X> (n) (2,*)Tj n , (3.5) 

n=0 

where the sum starts from 77 = as S(d,z), Eq. i|2.4[l . is regular at d — 2; the singularities in (d — 2) 
of the r.h.s. are then entirely due to the two last terms in the r.h.s. of Ea. l|3.3|l . with the double and 
simple poles in {d — 2) coming from the tadpoles entering in Eq. (|3.2() . That implies in particular that the 
Laurent-expansion of the l.h.s. must also have a double and a simple pole, as already anticipated. When 
expanding also the overall d-depending coefficient, one finds 

S(- 2 )(4,z) = -|, (3.6) 

S(- l \4,z) = ^ + 18), (3.7) 
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Note that the singular part of S(d, z) for d — > 4 is entirely determined by the above equations (and of 
course in agreement with previous results, see for instance Eq.s(53) of [1]). 

In the following of this paper we will show how to obtain S' (0) (2, z) and S (1 '>(2, z); 5 (0) (4, z) and S r(1) (4, z) 
can then be obtained by the previous formulae. 

4 The expansion around d = 2. 

It was shown in the previous Section, Eq. (|3.6l3.9|l that the coefficients S'^ n ' ) (4, z) of the expansion of S(d, z) 
around d = 4 can be expressed in terms of the coefficients S^ (2, z) of the expansion around d = 2, which 



is perhaps also simpler, as polar terms in (d — 2) are absent. In the following we will therefore restrict 
ourselves to the evaluation of the coefficients S^ n ' (2, z) of the expansion at d = 2, referring to Eo.s H3.6l3.9|) 
for obtaining the S^ (4, z). 

By inserting the expansion Eq. (|3.5(l into Ea. (|2.9|) and systematically expanding in (d — 2) also all the 
other d-depending terms, one obtains a system of chained equations, each equation corresponding to a 
given order in (d — 2). All the equations are of the form 

S (n) (2,z) = N ( - n \2,z) , (4.1) 

where n = 0, 1, .. is the order in the expansion in (d — 2). Note that the differential operator acting on 
S^ n '(2,z) is the same for any value of n, while the inhomogeneous terms do depend on n. The explicit 
expressions of the first two terms are 
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The system of equations Eq. (|4.1l4.2ll is chained, in the sense that the equation for the coefficient S^ n '(2, z) 
involves in general in the inhomogeneous term coefficients S^ ' (2, z) (and their first derivatives) of lower 
order k < n. Therefore the system must to be solved bottom up; as N^ a > (2, z) is explicitly known, one 
starts from n = and obtains S^°'(2, z), which then appears as a known term in the inhomogeneous part 
of the equation at n = 1 for S^ (2, z), and so on to higher values of n. 

The homogeneous part of Eq. (|4.1|l . which as already observed is the same for any value of n, reads 
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As it is a second order differential equation, it admits two linearly independent solutions, say ^i, (z), ^2, (2); 

^(*> = W^-W^, (4-4) 

is their Wronskian, according to Euler's method of the variation of the constants, the solutions of Eo. (|4.1|l 
are given by the integral representations 
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where $j , ^ are two integration constants. 

Up to this point Ea. H4.5[) is just formal: it will become a "substantial" , explicit analytic expression only 
when all the ingredients entering in it - the two ^i(z), their Wronskian and the integration constants ^ 
- will be actually evaluated. 

Although the Wronskian is defined in terms of the ^f%(z), it is known (Liouville's formula) that it can 
be immediately obtained (up to a multiplicative constant) from Eo. H4.3J) . An elementary calculation using 



the definition Eg. l|4.4|l and the value of the second derivatives of the ^i(z), as given by Eq. 14.3fl of which 
they are solutions, leads to the equation 

S^" — (j + ITT + iTi)"™' 

which gives at once 
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where the multiplicative constant has been fixed anticipating later results. 

Evaluating the two ^t(z) (and then the integration constants) is much more laborious. 

5 Building the solution of the homogeneous equation. 

Given a differential equation, it is immediate to obtain its solutions around any point zq in the form of a 
power series in (z — zo); indeed, when the solution is tentatively written as a power series, the differential 
equation can be easily recast as a recursive formula for the coefficients of the expansion. If the differential 
equation is of second order, as in Eq. (|4.3[) . the specification of the two initial conditions which identify the 
solution corresponds to the specification of the first two coefficients, needed for initializing the recursion. 

The expansion converges only up to the nearest singular point; in the case of Eq. (|4.3[) the singular points 
are the four points z = 0, — 1, —9 and z — oo. At each of the four singular points, one finds (see below) that 
Eo. (|4.3|l has a regular solution, say ip\(z), whose series expansion is straightforward, and a logarithmically 
divergent solution; the divergent solution can be written as the logarithm of the appropriate argument 
times the same function ij)x(z) which corresponds to the regular solution, plus another regular function, say 
ipi(z), whose series expansion coefficients are also recursively provided by the equation. 

While obtaining a local solution is immediate, building a solution in the whole range — oo < z < oo 
(i.e. working out the analytic continuation of a local solution) is a much more demanding task. At least in 
principle, one could consider unrelated pairs of linearly independent solutions evaluated as series expansions 
around some carefully chosen points, say Zj, each solution depending on a pair of constants (unrelated, to 
begin with), pick up one of the points, say zq, and start defining the two independent solutions at that point 
by fixing somehow the two constants at that point. One can then look at the expansion around the nearest 
point, say z\ (it is assumed that thanks to the careful choice of the points Zj the convergence regions of the 
expansions in zq and z\ overlap) and fix the two arbitrary constants of the expansion at z\ (i.e. obtain the 
analytic continuation from zq to zi) by imposing the equality of the solutions in the overlapping region. 
The procedure, which can be implemented without significant loss of precision also in a merely numerical 
way, can then be iterated to cover the whole range of interest of z. 

In order to cover with properly converging expansions also the singular points, which are at the boundary 
of the convergence regions of the expansions around nearby points, it is mandatory to consider also the 
expansions around the singular points themselves. Therefore, we will first work out the (still unrelated) 
solutions around each of the four singular points as series expansions; for joining them into a unique, 
analytically continued pair of solutions, rather than relying on expansions around auxiliary points (as in a 
numerical approach) we will look for "interpolating solutions" , valid within the intervals between successive 
pairs of singular points - namely the three intervals (0 > z > —1) , ( — 1 > z > —9) and (—9 > z > — oo) 
- and then will use the limiting values of the interpolating solution at the boundaries of each interval for 
relating the arbitrary constants of the solutions at the singular points, building in such a way the desired 
analytic continuation of the solutions in the whole range of the variable z. 

6 The solutions at the singular points. 

The presence of 1/z factors in the coefficients of Ea. (|4.3(l shows that z — is a singular point of the 
equation; when looking for a solution whose leading power in z is z a one finds the indicial equation a 2 = 0, 



indicating that the leading behaviors at z = are 1 and logz, so that around z = the two independent 
solutions can be written as 
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When -0^ (z),i/> 2 ( z ) are expanded as power series in z as 
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Eq. (|4.3|) gives for the coefficients V^ ( n ) the recursion relations valid for n > 1 
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determine uniquely the coefficients of the expansions. For definiteness, the first terms are 
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The radius of convergence of the series Eq. 1|6.2[) is 1 , as the nearby singularity of the solutions of Eq. (|4.3I) 
is z = — 1 or u = 1, see Ea. l|2.7|l . We were unable to solve the recurrence relations Ea. (|6.3(l in closed form, 
but by direct inspection one easily see that at large n the coefficients have a 1/n behaviour (which in turn 
implies that the radius of convergence is 1). For < z, i.e. > u (the Euclidean region), both solutions 
^l (z),^2 ( z )> Eq. (|6.1|l are real; as for spacelike u there are no singular points of Eq. (|4.3|) in the whole 
range > u > — oo, both functions can be continued along the whole positive z-axis < z < +oo, where 
they keep real values. For —1 < z < 0, if z = —{u + ie), with 1 > u > and e > 0, ^i (z) is still real, 
while 

In 2 = ln(— u — ie) = lnw — iir , 
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so that ^2 ( — u ~ * e ) develops an imaginary part equal to —i~Kip\ {—u 
Around z 



1, one finds similarly that two independent solutions of Eq. (|4.3(l are given by 
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where ip[ (z),ip 2 ( z ) are defined by the series expansions 
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with initial conditions 
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The recursion relations valid for n > 1 are 
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For definiteness, the first terms of the expansions are 
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As in the previous case, it was not possible to solve Eq.s l|6.9|) in closed form, but the coefficients are easily 
seen to behave as 1/n for large n, so that the radius of convergence of the series Ea. lj6.7fl is 1, the nearest 
singularity of the solutions of Eq.(^3J being at z = 0. Similarly to the previous case, ^\ (z) is real in the 
whole range of convergence of the expansion around z = —1, while ^\ (z) is real for z > — 1 and develops 
an imaginary part —iirifj^ (z) for z = — (u 4- ie), u > 1. 
Around z = — 9 one has similarly 
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where ip\ (z),ip 2 \ z ) are given by the series expansions 
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with the coefficients ij)\ „ determined by the initial conditions 
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and by the recursion relations valid for n > 1 



(») 



•01 






1 



72n 2 



72n 3 



(17n 2 -17n + 6)V>i 9 Li 
(17n - 12)4 9 Li 



-2(n-l)^_ 2 



(9) 



+n(17n 2 - Yin + 6)^2 „_i 



-n(n -1)V£L 2 



(6.14) 



For definitencss, the first terms of the expansions are 
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The expected radius of convergence of the series is 8 (the nearest singularity is at z = — 1 or u = 1), *fr\ (z) 
being real on the whole interval of convergence, while ^ 2 ( z ) develops an imaginary part —iwip\ (z) for 
z = — (u + ie), u > 9. 

Finally, for z — > 00 we put z = 1/y and find as independent solutions 
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with the coefficients V*i „ determined by the initial conditions 
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and by the recursion relations valid for n > 1 
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For definiteness, the first terms of the expansions are 

^{°°\z) = l-3y + 15y 2 
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The expected radius of convergence of the series in y = 1/z is 1/9, as the nearest singularity is at 
u = —z = 9. For y > 0, i.e. z > 0, the two series are real; for z = —(u + ie), oo > u > 9, hxy = — In it + i7r 
and $2 ( z ) develops an imaginary part —iir/u ipi (y). 

7 The interpolating solutions. 

Ea.s l4.H) . 14.2|) at n = give the equation valid for the n — 2 limiting value of S(d,z), Eq. H2.4fl : as the 
inhomogeneous term _/V(°)(2, z), Eq. l|4.2[l . is always real for real z, the imaginary part of S^(2, z) satisfies 
the homogeneous equation Eq. H4.3[) . As the imaginary part is present only above the physical threshold, 
i.e. when u = —z > 9, it is convenient to rewrite Eq. l|4.3[) in terms the function J(u) = ^(z) and of the 
variable u = — z, 



d 2 

lu 2 


'1 11 


d 

du 


1 1 1 


u u — 1 u — 9 


3u ' 4(u-l) 12(u-9)_ 



J(u) = 



(7.1) 



As already observed in in [9], the Cutkosky-Veltman rule gives for the imaginary part of <S'^ ^(2, z) (up to 
a multiplicative constant irrelevant here), the integral representation 



(v^-i) 2 



db 



y/Ri(u,b) 

where Ri(u, b) stands for the polynomial (of 4th order in b and 2nd order in u) 

R 4 (u, b) = 6(6 - 4)(6 - (V« - lf)(b - (Vu + I) 2 ) • 
For w > 9 the four roots in b of Ra{u, b) are ordered as 

o<4< (V^-i) 2 <(V^+i) 2 , 



(7.2) 



(7.3) 



(7.4) 



and the b integration in Eq. l|7.2[) runs between the two adjacent roots 4 and {\fu— l) 2 - We will now verify 
that Eq. I|7.2|l does indeed satisfy Eq. (|7.1|l by a procedure which will also suggest a broader set of similar 
solutions, also for other ranges of values of u. 

To that aim, let us introduce the auxiliary functions 



H(k,u) = 



dbb k 



y/R±{u,b) 



(7.5) 



/?, B are any pair of adjacent roots in b of R^iu^b), and k is an integer (if the integration boundaries do 
not include the root b = , k can take any positive or negative value; otherwise, k > for simplicity). One 
finds immediately 



db^ (b n y/Ri(u,b) 



= 



as the integration runs among two points at which the primitive of the integrand vanish; by working out 
explicitly the derivatives, moving always the square root y^R^u, b) to the denominator and using the 
definition Eq. l|7.5|l one finds the identity 



(k + 2)H(k + 3, u) - (2k + 3)(tt + 3)H(k + 2,u) + (k + l)(u + 3) 2 H(k + 1, u) 

-2{2k + l)(u-l) 2 H(k,u) = 



(7.6) 
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Another identity is obtained by writing 

B 



dby b(u-b + 3)-^IU{^b)J 

the derivative of the logarithm is 

-36+(m + 3) 

y/Ri(u,b) 

from which one finds 

3H(l,u)-(u + 3)H(0,u) = . (7.7) 

As a consequence of the identities Ea.s lj7.fil7.7jl . all the integrals H(k, it), Ea. (|7.5|l can be expressed in terms 
of only two "master integrals" , which can be taken to be 

H(0,u) ' "'' 



^R 4 {u 7 b) 

H(2,„) = / 7 H= ! . (7.8) 

J ^Ri(u,b) 

Also the u-derivatives of all the H(k, u) can be expressed in terms of the two functions. By using 
Ea.s H7.6l7.7|) write 

B 

Jdb^R^b) = (l + ^ u -~u 2 + lu^H(0,u)-(~ + ^H(2,u) 


B 

Jdbb^R^b) = (^ + ^u+^-u 2 + 7 -u 3 + ±u 4 + ±u^H(0,u) 

P 

15 3 9 o 1 ,\ rr ,„ . 

h -u u 2 u 3 H(2, u) . 

4 4 4 12 y V ' y 

One can then take the u-derivatives of both sides of the two equations; in the l.h.s. one can differentiate 
the integrands, and then express the result in terms of the two functions H(0, u), H(2, u) times polynomials 
in it; in the r.h.s. one finds both the two functions and their u-derivatives; solving for the derivatives gives 

±H M = Mu J ){u _ 9) [(3 + Uu-u 2 )H(0,u)^H(2,u)] 

—H(2,u) = — ^- \(9 + 228u + 30u 2 -12u 3 + u 4 )H(0,u) 

du ou(it — l)(it — 9) 

+ {-9-42u + 3u 2 )H(2,u)] . (7.9) 

With one further differentiation and some algebra one obtains also the second derivative 

^H(0,u) = 2u2{u - l ?{u - 9)2 [(-27 - 12, + 202^ - 36u' + u*)H(0, u) 

+ (27-60u + 9u 2 )H(2,u)] . (7.10) 

One can now substitute J(u) by H (0, u) in Ea. l|7.1[) : by using Ea.s l|7. 917. 10(1 for expressing the derivatives 
in terms of H(Q,u) and H(2,u), one finds that Eq. l|7.1|l is satisfied. But H(Q,u) Eq. (|7.8|l is equal to the 
integral Eq. (|7.2|) if j3 — 4 and B = {^u- l) 2 ; it is so verified that Eo. (|7.2|) satisfy Ea. l|7.8|l . as expected. 
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The verification that H(0, u) satisfies Eg. Q7.iJ [I remains valid for any choice of the two roots (/?, B) out 
of the four roots of Ri(u, b) listed in Eg. l|7.4|l . The linearly independent choices can be taken to be the 
three real integrals 

4 _ (v^-1) 2 ^ (v^+1) 2 

Jl{u) = 1 7 = m ' Mu) ^ / vTOT Mu) = I T^fer (7 - n) 

4 (v^-1) 2 

where the sign in front of Ri(u,b) within the square root has been adjusted to make all the integrals 
real (but an overall multiplicative factor i would be anyhow irrelevant when dealing with the solutions 
of a homogeneous equation). All the three functions Ji(u),i — 1,2,3 are solutions of Ea. H4.3f) . which has 
however only two linearly independent solutions; therefore they cannot be all independent. Indeed, consider 
the contour integral 

db - l , (7.12) 

y/Riiu, b) 

(oo) 

where the loop runs on the circle at infinity. For u > 9 the integrand has two cuts; one cut is from 
b = {^/u — l) 2 to b = (y/u + l) 2 and correspondingly its value is l/y/R^u, b) = l/{iy/—R±{u, b)) for 
b above the cut and — l/(ty— Ri(u, b)) for b below the cut, the other cut from 6 = to b = 4, and 
corresponding values — l/(iyJ—R±{u, b)) and l/(iW— R^u, b)) above and below the cut. The integral on 
the circle at infinity (where the integrand behaves as 1/6 2 ) vanishes; by shrinking the circle to two closed 
paths around the cuts, one obtains, apart an overall factor 2j, 

Ji(u) - J 3 («) = , (7.13) 

which reduces by one, as expected, the number of the independent solutions. 

The ordering Ea. l|7.4|l of the four roots in b of Ra{u, b) is valid, as already remarked, for u > 9; when 
u varies, the ordering of the roots varies as well, and the corresponding real interpolating solutions of the 
form of Eq. I|7.11(l are expressed by different choices of the integration interval in the variable b. 

8 The homogeneous solutions. 

We will now use the results of the previous Sections for building two independent solutions of Eq. (|4.3|l , say 
^i(z), ^2(2), properly continued in the whole range —00 < z < 00. We start from the the point z = 0; in 
the interval 1 > z > — 1 we define 

{tfi(z-ie) = *^ 0) (z-ie) , 
(8.1) 
V 2 (z-ie) = y { 2 0) (z-ie) , 

where the two functions ^[ (z), ^ (z) are defined by Ea. lj6.lfl and the equations following it. With that 
definition, the multiplicative constant of the Wronskian is fixed, and the Wronskian takes the value 

q 

W{z) 



z(z + l)(z + 9) ' 

already anticipated in Ea. lf4.fc)fl . 

Ea. l|8.1l) is meaningful within the convergence region of the series, i.e. 1 > z > — 1. For positive z the 
— it prescription is irrelevant, as the two functions ^i(z) are both real; for > z > — 1, i.e. < u < 1, 
^i(z) is still real, while the factor ln(^ — it) present in ^(z) develops an imaginary part —in, so that one 
has 

^ 1 {z-it) = ^f\z) , 
0>z>-l { (8.2) 

V 2 (z - it) = \n(-z)^\z) + ^ ( 2 °\z) - iwip[°\z) , 
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and, quite in general, for real z 

y i (z + ie) = (V i (z-ie))* . (8.3) 

In the z — > 0~ limit, in particular, one finds 

lim ty\{z — ie) = i , 

lim ^ 2 (z~ie) = ln(-z) - wr . (8.4) 

In the interval > z > — I, i.e. or < u < I, the four roots of R^u, b), Eo. (|7.3|) . are ordered as 

(0,(^-l) 2 ,(V^+l) 2 ,4) ; 
according to Sectional the interpolating solutions can be chosen as 

(v^-i) 2 

Jl (u) = / T^CT ' 

(v^+i) 2 



4 0ll) W = / , ^ • (8-5) 

The above functions are real, and regular in the region < u < 1; according to the discussion leading to 
Ea. (|7.I3[) . one has however the identity 

jf 1) (u) = J 1 ^( U ), (8.6) 

showing that only two of them can be linearly independent. In the z — > 0~ or u — ► + limit, the two roots 
(y/u — l) 2 and (y/u + l) 2 become equal, so that the integration in b of Eq.s (|8.5|) is elementary and gives 

lim j\ 0tl) (u) = -4= I - - hi M + In 3 

lim ji ' 1 ^) = -4= ■ (8.7) 

u-*o+ 2 V3 

The explicit calculation gives also 

lim jf' 1 '^) = -4= I --lnu + ln3 

in agreement, of course, with Eq.s l|8.6|l and (|8.7() . 

Given the different behaviours in the u — > + limit, J{ ' (m) and Jj ' (u) are linearly independent; in 
that same limit, one finds that their Wronskian 



lim 



jr\u)^\u)-±j^\u)jr\u) 



I 

6 ' 



shows the expected l/u singularity (the j\ ' '(u) are solutions of Eq. (|7.1|) . which is equivalent to Eo. (|4.3|) 
with z = —u and whose Wronskian Eq. (|4.6|l behaves as 1/z when z — ► 0, up to an overall multiplicative 
factor). 

The Vl'i(z), i = 1, 2, are linearly independent solutions of Eo. H4.3JI . so that the j\ ' (u), i = 1, 2, which 
are also solutions of the same equation (the equivalence of Eo. (|7.1[) and Eq. (|4.3fi has been repeatedly 
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recalled) can be expressed in terms of the ^i(z), and viceversa. By comparing the z — > behaviour of 
the &i(z) given by Eq.s l|8.4|l and the u — > + behaviour of the j] ' (u), Eq.s l8.7|l . one finds that in the 
interval > z > — 1, i.e. < u < 1 the solutions can be written as 

u = -z, > z > -1, < u < 1, 

tf 2 (*-*e) = -2V3J ( ,°' 1) (u) + — (21n3~nr)4°- 1) (u). (8.8) 

Similarly, in the u — > 1 _ limit, the two couples of roots (0, (\/u— l) 2 ) and ( (\/u+ 1) 2 , 4) become equal, 
so that one finds easily 

lim jf^) = lirn jf'^u) = \ , 

lim J 2 (0,1) (tt) = -~ln(l-u) + ~m2 . (8.9) 

u— >1— 4 4 

But for m — ► 1~ (or z — > — 1 + ) the two interpolating solutions j\ ' (u) can also be expressed in terms 
of the two ^f\ (z), Ea.s Hfci.6|) : matching the limiting values of Ea.s l8.9|l with the corresponding behaviours 
of Eq.s (|6.t)[) (the ^ (z) are real for < z < -1) gives 

J^\-z) = ^\z), 
0>z>-l { (8.10) 

jf 1) (-z) = lln2*«(z)-|vI/W(z). 

Substituting Ea.s H8.10|l into Ea.s l|8.8|l . we finally obtain 

*i(z - te) = ^ hi2*?\z - ie) - 9g*£\z - ie) , 

' ! ■'■' ' """" 2 *> ^2(z-ie) = ^(fln2ln3-2n-i9hi2)^\z-ie) (8 - H) 

+ ^ (-2 In 3 + in) ^ (z - le) , 

which gives the required analytic continuation of the solutions ^i(z), as defined by Eq. (|8.1|) . from the 
interval 1 > z > — 1 to the interval > z > —2 (or < u < 2), containing the singular point z = — 1 (or 
u = 1). 

The whole procedure can be repeated in the interval -9 > z > -1 or 1 < u < 9. The four roots of 
R±(u,b), Eq. (|7.3[) . are ordered as 

(o, (V^-i) 2 , 4, (V^ + i) 2 ), 

and the two interpolating solutions are given by 

(v^-i) 2 
J^\u) = ' '"' 







y/-Ri(u,b) 



J^\u) = f — d ±= , (8.12) 

J wRi(u,b) 

(v^-i) 2 



while for the third solution, according to Eq. l|7.13[l . one has 

(%Ar+i) 2 

J y f -Ri{U,b) 

4 



14 



In the u — » 1 + limit (as in the u —> 1 limit) the pairs of roots (0, (\A* — I) 2 ) an d ( {y/u+ 1) 2 ,4) 
become equal and one evaluates easily 

■n 
4 ' 

r(l,9), 



Urn Jf' 9) (w) 



lim J^(u) = -^ln( W -l) + jln2 



(8.13) 



As the J> ' J (u) are solutions of Eg. (|7.1fl . equivalent to Eg. lf4.3fl . in the interval 1 < u < 9 or — 1 > z > —9, 
they can be written as linear combinations of the &i(z— ie); the comparison of Eg.s lf8.13f l with Eg.s (f6.bl6.1Uf) 
in the z — ► — 1 _ limit gives 



' *i(z-ie) = ^(4 1 ' 9 \u) + 3ij[ 1 > 9 \u)) , 



-1 > z > -9 < 



v^S 



(1,9), 



tf 2 (z - ie) = ^ ( (tt + 6iln3)J} x, "'(«) 

+ (21n3-i7r)J^' 9) (u) 



(8.14) 



We look now at u — > 9 ; in that limit, (y 7 ?! — l) 2 — > 4, and one obtains 



lim Ji (1,9) (u) 

u->9~ 



4\/3 



(-m(9-u) + 31n2 + 21n3) 



lim J, ' (u) = — = 
u-9- 4\/3 



.15) 



In close analogy to what already done for the interval < u < 1, in the interval 1 < u < 9 the two interpo- 
lating solutions J> ' (u) can be expressed in terms of the ^\ (z) Eo.s (|6.11|l : matching the coefficients of 
the linear combination by using the limiting values for u — ► 9 — , one obtains 



-1 > z > -9 



J i (") = 173 ((31n2 + 2]n3)*H*)-*W(*) 



(8.16) 



When substituting Eg.s lj8.16fl in Ea.s f8.14fl one finally obtains 



*!(«) = i [tt + i 3(3 In 2 + 2 In 3)] ^ a; (z) 



(9), 



3 , T ,(9) 



-1 > z > -17 { 



*&*n*)> 



* 2 («)= [|ln2 + ln3+^(181n21n3 + 121n 2 3-7r 2 )]*i 9) (z) 



.17) 



-^(^ + z61n3)^ 9) (z) 



>(»)/ 



The radius of convergence of the expansion around z = —9, i.e. u = 9, of the \£^ (z) is 8, as the nearest 
singularity is at u = 1 (or z = —1), hence the range of validity — 1 > z > — 17orl<u<17of the above 
formula. 

The last interval to consider is 9 < u < oo; the four roots of Ra(u, b) are ordered as 

o,4,(V^-i) 2 ,(V^ + i) 2 , 

and two interpolating solutions corresponding to integrating in b between two adjacent roots are 



j[ 9 '°°\u) 



jJ 9 ->(«) 



db 



y/-Ik(u,b) ' 



(v^-1) 2 



db 



y/R 4 (u,b) 



(8.18) 
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while the third solution satisfies 

(xAT+l) 2 

J y/-Ri(u,b) 

(v^-l) 2 

Proceedings as int the previous case, from the limiting values for u — > 9 + 

lim jf ,oc) (u) = — =(-ln(u-9)+31n2 + 21n3) 

u->9+ 4v3 



lim + 4 9 ' oo) (u) = -^=. (8.19) 

«^9+ 4V3 



one obtains 

' *i(* - ie) - ^ (-2jf°°>(u) + 3zjf ' oo) («)) 



-9 > z > -oo < 



* 2 (z - ie) = *fi ( (tt + 6iln3)Jf ,uu » 
-41n3J 2 (9>oo) (u)) . 



(8.20) 



The limiting values for u — > oo are 

lim </} '°° (u) = — 



u 
hm ji 9,oo) (w) = ^-lnu, (8.21) 

implying 

f 4 9 -° C \u) = ~n^r\z- i e), 
-9 > 2 > -oo < (8.22) 

[ J^'°°)(tt) = -li^[°°\z - ie) + ^t\z - ie) . 

Note that according to the definitions of the ^f\°° (z), Eq. (|6.1t)[l . the r.h.s. of the above equations are real, 
as expected (according to the definition Eq. l|8.18[l . the j\ ' (it) are indeed real). 
By substituting Eq.s (|8.22|l into Eq.s (|8.20|l one finally obtains, 

%l(z-K)=-$&* ( £°\z-ie), 

-9 > z > -oo { (8.23) 

* 2 (z - ie) = -Vin^iz ~ ie) - 6V3^r^Hz - ie) . 

The formula holds also in the range oo > z > 9 (where the ie in the arguments can be ignored). 

9 Identities under transformations of the argument. 

Consider the transformation 

9 9 

y=~, z=-. (9.1) 

z y 

It maps the singular points of Eo. (|4.3() . z = 0, — 1,— 9,oo into themselves, more exactly — > oo — » 0, 

(—1) — > (—9) — ► (— 1). If z is real and negative, y is also real and negative, and the ie prescription is 

implemented as 

9 9 9 9 , . 

= — h ie = y + ie , = ie = z — ie . (9-2) 

2 — ie z y + ie y 

As it is easy to verify, if 'I'(z) is a solution of Eo. H4.3fl . 

i. (|) (9.3) 
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is also a solution of the same equation, and can be therefore expressed in terms of the solutions ^i(z) and 
$2(2) defined in the previous Section. The argument applies, of course, to the ^i(z) themselves, for which 
we can therefore write 

*< ( I -ie) = \z V AijVjiz + ie) , (9.4) 

KZ J ^ 3=1,2 

where the coefficients Aij are the elements of a 2 x 2 matrix A and the overall factor 1/3 has been introduced 
for convenience. One can obtain the four coefficients Ay by imposing for instance Eq. (|9.4l) to be valid for 
z = — 1 + r), y = —9 — 77, with 77 — > + , were the limiting values of the ^Si(z) are known from Eci.s H6.6l6.llf) . 
The result is 

A ^> ( 2hl3 ^ ^ «XK\ 

A= — \A\n^-\^ -21n3 ) ' (9 ' 5) 

One can check that with the above values of the Aij Ea. l|9.4|l holds also in the z — > limit, when the 
limiting values of the ^i(z) are given by Eq.s (|6.1l6.16|) . 
Ea. l9.4|) can also be written, letting z — » 9/z, as 

^(*-ie) = 3~ J2 4y*i(~ + 

Z 3=1,2 \ z 

as the matrix A is real, by taking the complex conjugate and recalling Eo. (|8.3|) one has 



*<(« + ie) = 3 - J2 AijVj ( - - *c) , (9.6) 



z ' — ' \ z 

3=1,2 v 

When chaining Eq. (|9.6(l and Ea. H9.4|l one finds the condition 

A-A=l, (9.7) 

which is of course satisfied by the explicit values given in Eq. (|9.5|) . Eq. (|9.7|l shows also that two of the four 
coefficients Ay are fixed once the other two are given. 

The transformation Eq. l|9.1|l maps the points z — —3, 3 into themselves. Correspondingly, Eo.s H9.4l9.5(l 
give 

*a(3) = ^(2V31n3-7r) *i(3) , 
v3 

i?e* 2 (-3) = -^(2V31n3 + 7r) fle#i(-3) , 
V3 

Jm* 2 (-3) = -p(2-s/31n3-7r) 7m*i(-3) . (9.8) 

v3 

In terms of the interpolating solutions, Ea.s H9.4[l become: in the interval < u < 1, with 00 > 9/u > 9, 

0< u < 1 

,(9,oo) /9\ _ V3 ( o,i) 



Jt '[^ = ^r»; (9-9) 



in the interval 1 < u < 9, with 9 > 9/u > 1, 

1 < u <9 

.(1,9) f 9 \ _ ^3,(1,9), 



J * uj = «x Ji (ti) ' (9 - 10) 
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and, in the interval 9 < u < oo, with 1 > 9/u > 0, 

9 < u < OO 

A™ (I) = ufj^Hu), 

4° A) (l) = ufj^\u). (9.11) 

Ea.s H9.11|) are of course equivalent to Eg.s (|9.9[) : similarly, the second of Eg.s l|9.10|l is equivalent to the first. 
At u — 3, in particular, they reduce to 

J^ 9) (3) = ViJ { ^ 9 \i) . (9.12) 

We consider next the transformation 

y = — — T - z = — — T • (9-13) 

z + 1 y + 1 

which also maps into themselves the singular points of Eo. (|4.3l) . more exactly it maps — ► (—9) — > 0, and 
(— 1) — ► oo — ► (— 1). For real y, z the transformation implies for the ie prescription 



z — ie + 9 
z — ie + 1 


z + 9 
z + 1 


- ie = y - ie , 


y — ie + 9 


2/ + 9 


- ie = z — ie . 



(9-14) 
y — ie + 1 y + 1 



One can verify that, if \P(z) is a solution of Eg. (14. 31) . 



1 .tf-i±? , | (9.15) 



z+1 V z + 1 

is also a solution of the same equation, and therefore a combination of ^i{z) and \&2(z)- When applied to 
the ^i(z) themselves the argument gives 

^ Z ' .3=1,2 

where the coefficients i?^ are the elements of a 2 x 2 matrix B and the overall factor v2/4 has been 
introduced for convenience. The Bij can be fixed, for instance, by imposing the validity of Eo. (|9.16ll for 
z — > 0; the result reads 

V2/ 7r + i61n3 -»3 \ , . 

2tt V 47rln3 + i(121n 2 3-7r 2 ) -7r-?:61n3y' l ' 

One can check that with the above values of the B^ Eq. (|9.16|) holds also in the z — ► — 1 limit. 
Substituting z by —(2 + l)/(z + 9) Eq. l|9.16[l can also be written as 

2a/2 v- / ^ + 9 \ , 

*i(* - »e) = -j^t E %*j V ITT " l 7 ; (9 ' 18) 

j = r,2 

by chaining Eq. I|9.16|) and Eq. (|9.18|l one finds 

B-B = -l, (9.19) 

satisfied of course by the explicit values Eq. (|9.17|l , showing again that two of the four coefficients B^ are 
fixed once the other two are given. 
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The transformation Eo. (|9.13|) maps z = 3 into z — —3. Correspondingly, from Ea.s (|9.16l9.17|) and 
Ea. (|9.8|l one obtains 

Re^ 1 (-3) = *x(3) , 
/m*i(-3) = -V3*i(3) , 

i?etf 2 (-3) = -(V37r + 61n3)*i(3) , 

7m* 2 (-3) = (7r-2\/31n3)*i(3) . (9.20) 

In terms of the interpolating solutions, in the interval < u < 1, with 9 < (9 — u)/(l — u) < oo, 
Ea.s l9.16l) become 

< u < 1 



r (9,oo) { 9 ~ u \ l ~ u t(0,1)/ x 



•# ,00) (l^) = ^ 0,1) («). (9-2D 

The range 1 < u < 9 is mapped by the transformation into the (spacelike) interval — oo < — (9— u)/(u— 1) < 
0, where we did not consider the interpolating solutions; the range 9 < u < oo, finally, is mapped into 
< (u - 9)/(u - 1) < 1, where Eq.s (|9.16[) give 

9 < u < oo 



r (0,l) ( u -$ \ = u ~ l T (9,' 

u - 1 / 2 



J^r—i) = ^J^iu). (9.22) 



The above relations are of course equivalent to Eq.s (|9.21|l . 

The transformations Eq.s l|9.1l9.13[) can be combined into a third transformation 

z+1 y+1 , 

V = -9 -±r , z = -9 2— , 9.23 

z +9 y + 9 

which maps the singular points of Eq. l|4.3|l into themselves, — ► — 1 — ► 0, and —9 — > oo — > —9. For real 
y, z it implies for the ze prescription 

z — ie + 1 z + 1 

- 9 7^7T9 = " 9 7T9 + 4£ = y + i£ ' 

_ 9 ^i±i = _ 9 y±i + . e = ^ + . e (924) 

y — ie + 9 y + 9 



Again, if ^(z) is a solution of Eq. (|4.3ll . 



1 / z + T. 

n* "9 9.25 

z + 9 V z + 9 ' 



is also a solution of the same equation, as it is easy to verify, and therefore can be expressed in terms of 
^i(z) and * 2 (z) as 

** (~ 9 Jrl + i£ ) = if (z + 9) S c ^ z - ?£ ) ' ( 9 - 26 ) 

where C is a 2 x 2 matrix, and the overall constant \/2/12 has been introduced for convenience. 
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By writing Ea. H9.6fl with z replaced by — 9(z + l)/(z + 9) and then using Eq. lf9.16f) one finds 

t ( n Z + l ■ \ lz + 9 v^ , _ / z + 9 

V ^ + 9 / 3z + l £f v J \ z + \ 

= ~^§( Z + V) E E AiiBjkVkiz-U). (9.27) 

3=1,2 fc=l,2 



V2 V3 / 21n3-i7r -1 



By comparison with Eg. (|9. 26(1 . 

C = -A-B = ±^-^[ "To i ' • n1 " • ) ■ ( 9 - 28 ) 

2?r ^lnS + ^Tr 2 -21n3-i7ry v ; 

Exchanging (z — ie) with — 9(z + l)/(z + 9) + ie, Eg. (|9. 261 can also be written as 

j — lj-^ 

Note the appearance of C* in the r.h.s of the previous formula; indeed, the sign of ie in the arguments has 
changed with respect to Ea. H9.26|) . and according to Eq. H8.3fl that amounts to take the complex conjugate. 
By chaining Eo. (|9.29fl and Ea. H9.26fl one gets the identity 

C* ■ C = 1 , (9.30) 

which is of course satisfied by the explicit values of the coefficients CV,- given by Eq. (f9.28fl . In terms of the 
interpolating solutions, in the interval < u < 1, with 1 > 9(1 — u)/(9 — u) > 0, Ea.s lf9.26|) become 

< u < 1 

r(0,l) ( a ±-JL\ - Ao_.,w(°>l) 






J * [ 9 9^) = l- (9_u)Jl ( " ); (9 - 31) 

the interval 1 < u < 9 is mapped into the spacelike range > 9(1— u)/ (9— u) > — oo, where the interpolating 
functions were not considered, while for 9 < u < oo, corresponding to oo > 9(1 — u)/(9 — u) > 9, one finds 
the relations 

9 < u < oo 

^~ ) (^) - i|(«-9)^ ,0O) («); 

r(9,oo) /ott-l^ _ V3,.,_ Q w(9,oo), 



J 2 ^^J = lf(«-9^i '(«)> ( 9 - 32 ) 

which arc of course equivalent to Eq.s (|9.31fl . 

10 The solution S'^(2, z) at zeroth order in [d— 2). 

We can now reconsider Eq. (|4.5fl for n = 0, i.e. at zeroth order in the expansion in (d— 2). With the explicit 
value of 7V(°)(2, z) given by the first of Eo.s l4.2fl . Eo. lf4.5fi becomes 

SV>(2,z) = *i(z) U® - ± J dw * 2 (w) 

+ *„(*) l*V> + ±jdw*i(w)) , (10.1) 
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where the functions ^i(z) are known from the previous Sections, while the two integration constants are 
still to be fixed. To fix them, we will exploit the information that the solution, corresponding to the loop 
integral Eq.(refeq:defMI), is real as far asw = — z is below the physical threshold u = 9, i.e. in the whole 
range — oo<w<9oroo>z> — 9, with the reality condition valid for any d, in particular for the zeroth 
order term of the expansion around d = 2. 

In the range > z > —1, i.e. < u < 1, using Eq.s (|8.8|l for expressing the *&i(z) in terms of the 
Jl ' (u), the reality condition gives 

Im S (0) (2,z) = -V34 0A) (u) ¥° } , (10.2) 

which implies 

* 2 0) = . (10.3) 

In the range —l>z> —9, i.e. 1 < u < 9, using Eci.s H8.8|) and Ea.s l8.14|) one obtains 

Im S<- \2,z) = I I VM 0) - \ fdv Jf'%) ) j[ x '»\u) ; (10.4) 

the (qualitative) condition that S^(2,z) is real therefore gives the (quantitative) result 

^^jdoJi^iv), (10.5) 



or, on account of Ea. lB.5|) of the Appendix, 

v/3 



•f-^Cs)- <">*> 



where Ci2(</>) is the Clausen function of weight 2. 

Ea.s (|10.3ll0.()|l determine completely the solution Eq. (|10.1|) . which from now on can be considered 
as known in closed analytic form. Indeed, expressing again, when needed, the ^i(z) in terms of the 
interpolating solutions Ji (u) , and on account of the results of the Appendix, we can evaluate in particular 
the behaviour of S^°\2, z) at the singular points of Ea. (|4.1|l at n = 0. We find: 



• at z = 0, 



l 



* (0) (2,0) = *< 0) = ^Jdv ^(v) = f Cl 2 (|) , (10.7) 



in agreement with a result of [10]; 
• on the mass shell — z = u = 1 



o 
at the threshold (z + 9) -> 0+ or -z = u -> 9~ 

'z + 9 



± 

S(°)(2,-l) = i|AJr»( B ),i, 2 ; (10.8) 



S<°>(2,z) ► -^ 

v ; ( 2 +9)^o+ 48 



7T 111 



+ 5C1 2 (- 



..., , , .. y . tl +0(z + 9) ; (10.9) 

above threshold, with z — — (u + ie) and u > 9, S , '°'(2, z) develops the imaginary part 

ImSW(2,-u-ie) = ^ ./^(u) , (10.10) 

in agreement with the comments accompanying the introduction of Eq. (|7.2|) and Eq. (|8.18|l : 
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• finally, when z — » +00 (spacelike region), the behaviour of S^ (2, z) is 

^°)(2,z) — -— -|- \n 2 z + o(\\ . (10.11) 

z^+ixj 16z \ z J 

The above results show explicitly that S^\2, z) is regular at z = and z = — 1, as already recalled. It 
can be interesting to observe that by performing the change of variable v — > 9/u, and then using the first 
of Eq.s (|9~TTjl . Ea. ((Tir7JI becomes 

OO 

5(0)(2,0) = i|^jf oo) («); (10.12) 

9 

similarly, with the change v —> (v — 9)/(v — 1) and the second of Ea.s l9.22|) . Ea. l|10.8|l reads 

00 
S<°>(2, -1) = - / -^r Ji 9,oo) («) • (10.13) 



47 v-1 2 

9 

Quite in general, one can write for 5^°- ) (2, z) the dispersion relation 

00 00 

S^(2,z) = ± l^-ImS«\2,-v-ie) = \ f J?-jl*>~\v) , (10.14) 

n J v + z 4 7 V + z 

9 9 

where Eq. (|10.10fl has been used; Eq. (|10.12|l and Ea. H10.13() can then be seen as the dispersion relation 
evaluated at z = and z = — 1. 

We can also easily work out the expansions of S^°'(2,z) around any of the singular points points of 
Ea. (|4.1(l at n = 0, namely z — (0, — 1, —9, 00), with the above results providing the initial values and using 
either the known expansions of the ^i(z) given in Section |S| or, better, the differential equation Ea. (l4.1|) 
for obtaining the coefficients of the expansions up to any required order. We find: 

• around z — 0, according to Ea. (|10.7J) . the expansion can be written as 

5^°)(2,z) = *( 0) + J2 4 0)zfc ; ( 10 - 15 ) 

fc=ljOO 

substituting in Eq. jO]| with n = 0, one finds 

S (0;0) M = ^(0)^(0) {z) 

+ —z-—z 2 + ...., (10.16) 

24 864 ' V ' 

where ip\ {£) is the same as in Ea. (|6.2|l ; 

• around z = — 1, according to Eq. (|10.8f) . the expansion can be written as 

S«W(2,,) = £ + £ -L 1) (* + l)*. (10-17) 



64 

k=i,. 



and the differential equation gives 



S (0;1) M = ^\z) 



&{* + !)- ±{z + !)> + .... (10-18) 



<,(i), 



where ip\ (z) is the same as in Eo. (16.711 ; 
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around z — —9, according to Eq. (|10.9fl . the expansion can be written as 

S^>(2,z)=ln(z + 9) £ 4 9) (- + 9) fe + E 4 9) (- + 9) A 

k=0,oo 



k— O.oo 



with 



.(9) 



48"' 



the differential equation gives 

S(° ;9 )(2, 



.(9) 



,(9) 



V3 
48 



-7rm(72)+5Cl 2 (| 



,,(9) 



Hz + 9)^r(z) + ^'(z) 



(9), 



+ 1^2^ + 9 ) + 69T2( Z + 9 ) 2 + -- 
where the tf)\ (z) are the same as in Eq. (|6.f2|l ; 
for large z, finally, according to Eq. (|10.11|) the expansion can be written as 



S^' oo) (2 1 z) = - 



^ E -t O0) i+^ E *ir°i+ E 4°< 



/c— 0,oo 



with 



,(°°) 



fe=i, 
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fe=i,i 



Imposing the validity of the differential equation, the expansion can be written as 



S (o;«0 M = ^_ 

loz 



In 2 z ^°°> f-^j 



— 2 In z %j), 



(oo) 



2 1 

~z ~ 2~? 



where the ?/^ (z) are the same as in Eq. Htj.l7|) . 



(10.19) 



(10.20) 



(10.21) 



(10.22) 



(10.23) 



The correspoding expansions for dS^ (2, z)/dz is immediately obtained by differentiation. From Eq. (|2.10|l 
one can then obtain the z-expansions for S\ (2, z), the zeroth order term in (d — 2) of the master integral 
Si(d,z) Ea.i2~5l). 

The above z-cxpansions can be used as the starting building blocks for implementing a computer routine 
for the fast and precise numerical evaluation of S^(2, z). 



11 The solution S^\^z) at zeroth order in [d — 4). 

From the results of the previous Section and from Eo. (|3.8|) one immediately obtains the relevant z- 
expansions for S , ^ ^(4, z), the zeroth order term in (d — 4) of S(d, z) Eq. l|2.4|l : 



• around z = 



S(0;0) (4 } = ^(0)^3 1 J_ 

v ' ' l \2 3 27 
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11 



32 128 1728 



z z + 



(11.1) 



v(0) 



with ^\ ' given by Eq.dHl; 
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around z 



S^(4,z) = Zj! (-!(* + l)*-jj(»+ 1)3 + ....) 



128 + 128^ + 1 ) + 64 (2 + 1) + 384 (Z + 1) + 



(11.2) 



• around z 






S-^ +9 )-^ +9 ) 2+ 



10368 



(z + 9) 3 + 



with 4 9) and t^ ] as in Eg.l |10.20f> : 
• for large and positive z (spacelike u), finally 



5 (0;oo) (4z) = ln 2 / 3 3 



32 16z 16z 2 



, ( 1 9 3 

In z — z H 1 + 

V32 32z 8z 2 
13 15 3 29 

128 Z ~ 32 ~ 64i + 64^ 



(11.3) 



(11.4) 



The above expansions can be used as the starting building blocks for the fast and precise numerical evalu- 
ation of 5^(4, z). 



12 The solution 5' (1) (2, z) at first order in (d-2). 

Once S^(2, z), is known, N^(2, z), given by the second of Eq.s (|4.2() . is also known, so that we can work 
out Ea. l|4.1|l and its solution Ea. l4.5|l for n = 1, i.e. at first order in the expansion in {d — 2). After some 
algebra, it takes the relatively simple expression 



S«(2,z) = * x (z) 



* 



(i) 



v(0) 



, , *i u; ( - In 3 - - lnz + - ln(z + 1) + - ln(z + 9) 



1 /". L ! 2 2 

— I dw[A 1 1 

96 / V w w + 1 w + 9 



* 2 (w) 



+ ^ 2 (z) 



* 



(i) , 1 i r,(o) 



-* 



1 f f 1 2 2 

— dw \ A 1 1 

96 J \ w w + l w + 9 
o 



*i(tu) 



(12.1) 



where ty\ ' is given by Ea. pU.6|l (according to Ea. (|10.3[l ^> 2 Q) an d the two integration constants $f\ 
are still to be determined. Imposing as in the previous Section the reality conditions for oo > z > —9, we 
find 

(12.2) 



*« = 



24 



and 



*« 



73 
24 

12 



rfc 



|j| 0,1) (v) + ( 2 + In 3 + i Inv - ln(l - ») - ln(9 - v) ) j[ ' h ( r) 



CI; 



(i) 



3 /27T 

i LS3 (t 



1 



ln3Cl 2 



(D-a- 1 



(12.3) 



where use is made of the definite integrals of the Appendix. 

S'( 1 ' ) (2, z) is now fully determined in closed analytic form. Expressing again the ^i(z) in terms of the 
interpolating solutions Ji(u) and on account of the results of the Appendix, we obtain in particular the 
following behaviours: 



• at z = 



S<.> ( 2,0) = <» = ^| 



*©-^!^«5)4-' 



(12.4) 



in agreement with [10] if allowance is made for the difference in our integration measure Eq. l|2.1|l . 
which generates the extra term Cl(7r/3); 



• on the mass shell — z = u = 1 

S (1) (2,-l) 



1(^ + 3^1112-^(3) 



(12.5) 



at the threshold [z + 9) — > + or — z = u — ► 9 
SW(2,z) 



(z+9)^0+ 



- ^ lnfi±^ 
96 V 72 



V3 



64 



Ls 3 



2tt 



ln(z + 9) + 31n2-ln3 + 2 
l(2 + 61n2 + ln3)Cl 2 (^ 



finally, when z — > +oo (spacelike region), the behaviour of S^(2, z) is 



S (1) (2,z) 



1 



z^+oo 32z 



3 In 3 z + 6 In 2 z - 7T 2 In z + 18C(3) 



O 



(12.6) 



(12.7) 



The corresponding complete expansions in z for S^'(2, z), dS < - 1 '(2, z)/dz, S ,< - 1 '(4, z) etc. can be obtained 
as in the case of S^°'(2, z); we omit them for the sake of brevity. Let us just observe that from Eq. (|3.9|) 
and Ea. i|12.5fl one gets 



Go 



G. r ) 



S<"(4,-D = |^(2,-l) + ^ = |j + ^,«*i S (l,4,-l) 



(12.. 



in agreement with the known result [11]. 

13 Conclusions. 

We have considered the system of first order differential equations in the external momentum transfer 
satisfied by the Master Integrals of the two loop sunrise graph in the equal mass limit and in the usual 
continuous d-dimensional regularisation. The system is equivalent to a single second order equation for 
the scalar MI; after recalling the relation between the expansions in (d — 4) and (d — 2), the second order 
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equation is expanded around d = 2 and solved by means of the variation of the constants formula by Euler. 
That requires the knowledge of the solutions of the homogeneous equation, which is obtained by working 
out first the series expansions of the differential equation at all its singular points, then joining smoothly 
those expansions by means of a set of "interpolating solutions" , whose explicit form (one-dimensional 
definite integrals of a suitable integrand) is suggested by the Cutkosky-Veltman rule for the imaginary part 
of the MI. The interpolating solutions, whose relation to the complete elliptic integrals of the first kind 
is discussed in an Appendix, are found to transform on themselves for a particular set of transformations 
of the arguments. Once the full knowledge of the homogeneous equation is established, Euler's formula 
provides in closed analytic form also the limiting values of the solutions at all the singular points of the 
equation. They are worked out explicitly for the zeroth order in (d— 2), the related zeroth order in (d— 4) 
and the first order in (d — 2). 

The analytic behaviour at the singular points and the related expansions can be used as the starting 
building blocks for the implementation of a fast and precise numerical routine for the evaluation of the 
equal-mass sunrise scalar integral for arbitrary values of the momentum transfer. The extension of the 
approach to other two-loop self-mass amplitudes with a single mass scale seems easy to obtain; in that way 
one could complete the analytic and numerical evaluation of the two-loop electron self-mass in QED. While 
such a result is somewhat academical, much more interesting, even if correspondingly difficult (but not out 
of reach) will be the extension of the approach to the two-loop self-mass integrals in the general mass case, 
as required by the precision tests of the current Standard Model (QCD & EW interactions) of elementary 
particles. 
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A Relations with the Complete Elliptic Integrals. 

The integral Eq. (|7.2l) . which is also the phase space of three particles of equal mass, can be written as 
a Complete Elliptic Integral of the First Kind [9], (the first explicit reference to phase space integrals 
and elliptic integrals was perhaps made in [13]); the same is true of course for all the other integrals 
corresponding to the interpolating solutions. 

Quite in general, given a polynomial of fourth order in b, written as 

R 4 (b) = (b- 6i)(6 - b 2 ){b - b 3 )(b - b 4 ) , (A.l) 

where the four roots are ordered as 



and the integral 



&i < b 2 < b 3 < b 4 , (A. 2) 



J(6 ll 6 2 )= / -r^=, (A.3) 



bl 



performing the standard change of variable [14] 



one obtains 



a (b4-b 2 )(b-b 1 ) b 4 (6 2 -b 1 )x 2 + & 1 (6 4 -& 2 ) ,. . 

X (b2-h)(h-b)> (h - h)x2 + (& 4 - b 2 ) ' [ ' 



bi 



2<> 



where 

(64 - 62) (63 - 61) 

and K(m) is the Complete Elliptic Integral of the First Kind, 

1 

K(m) = f dX = . (A.7) 

J J(l-x 2 )(l-mx 2 ) 


Similarly, one finds 

,{hM -]iim= 7^m^T) K{i - n) ' (A ' 8) 

where m is the same as in Eq. (|A.6|) . and finally 

I{hM -l7vr 7^m^o K{mh (A ' 9) 

so that, quite in general 

7(61,62) =1(63,64) , (A.10) 

in agreement with Eg. s (|7. 1117. 13(1 . 

By using the above results, we can express the interpolating solutions introduced in Section |H1 in terms 
of complete elliptic integrals. We find: 

• for the interpolating solutions in the interval < u < 1, defined by Ea.s l|8.5|l . 

j| 0,1) (u) = 1 K(a(u)) , 

v/(i + V^) 3 (3-V^) } 

4 0,1) (u) = 1 K(l-a(u)), 

v/(i + V^) 3 (3-V^) 

(l-» 3 (3 + 0I) 

a(u) = n , ^ 3 , Q 7=7 ; (A.ll) 

(1 + V") (3 - V - ") 

• for the interpolating solutions in the interval 1 < u < 9, defined by Ea.s H8.12[) . 



j! "' s,( » ) - 7h K(m ■ 

J«'' 8| («) = -J=A"(i -&(„)), 

Vlovu 

_ bfrjtea (A12) 

lOy/U 

finally, for the interpolating solutions in the interval 9 < u < 00, defined by Eq.s (|8.18|l . 

^l 9,00) («) = X : =K(C(U)) , 

■^M = /rA - ' ^(1 - c(«)) , 

c(u) = ^ ■ (A.13) 

(V - "- !) (3 + V u ) 
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The identities under transformations of the argument, written in Section[!5]for the interpolating solutions, 
can be rewritten in terms of the elliptic integrals. We obtain, in terms of the arguments a(u) , b(u) , c(u) 
defined in Eq.s (|A.lllA.12IA.13|) : 



• from Ea.sl(0| or Eq.s l^TTT]) 

with t = 9/u; 

• from Eq.sljnUni) 

again with t = 9/u; 

• from Ea.s ij9~22|) or Eq.s (|^T|l 



with v = (u-9)/(u- 1); 

• from Ea.s l|9~3T|) 

with w = 9(u- l)/(tt-9); 

• finally from Eq.s(| 



K(l-a(u)) K(c(t)) 

K(a(u)) K(l-c(t)) ' l ' 



K(l- b{ u)) =3 K {m ) 



K(b(u)) K(l-b{t)) ' 



K(l-a(u)) =2 K(l-c(v)) ^ (Ai6) 



K(a(u)) K(c(v)) 



K(l-a{u)) K(a(w)) 

K(a(u)) K(l-a(w)) ' l ' 



K(l-c(u)) _ 3 K(c{w)) 



K(c(u)) 2K(l-c(w)) 



(A.18) 



again with w = 9(u — l)/(w — 9). 



Eq. ljA.14p is based on the algebraic transformation a(u) — > 1 — c(9/u); when the parameter u is in 
the range < u < 1 all the arguments are real and in the range (0, 1), and all the elliptic functions take 
real and positive values (the equation remains of course valid, by analytic continuation, for any value of 
u). Similarly, all the arguments are in the range (0,1) in Eo.s (|A.16|l and l|A.17(l also for < u < 1, in 
Eo. (|A.15|l for 1 < u < 9 and for 9 < u < oo in Ea. (|A.18(l . but the relations can be continued to arbitrary 
values of u. The algebraic relation derived from Ea. (|A.14(l by eliminating the parameter u between a(u) 
and 1 — c(9/u) is a modular equation of degree 3 [15]; similarly, Eq. (|A.15|) gives a modular equation also 
of degree 3 between b(u) and 1 — b(9/u), Eq. (|A.16|) a modular equation of degree 2 between a(u) and 
c{{u — 9)/(u— 1)), Ea. (|A.17[l an equation of degree 6 between a(u) and 1 — a(9(l — w)/(9 — u)) and finally 
Eo. (|A.18|) a modular equation of degree 3/2 between c(u) and 1 — c(9(l — u)/(9 — u)). 

B Definite Integrals. 

In this appendix we discuss the evaluation of the definite integrals used in the previous Sections. As a first 
case, we consider 

i i (i-V^) 2 

duJi°'%)=fdu / - = f = , (B.l) 

J J \/ —R±{u,b) 



2,s 



where j\ ' ' (u) is defined by Ea. lj8.5fl and R±{u, b) by Ea. f7.3ft . One can invert the order of integration, 
obtaining 



duJ^ 1] {u) 



db 



(i~Vb) 2 



du 



n VK^-b) I ^ {{ i-Vbf-u){{i + VbY-u) 



1 



db 



In b 
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where the u-integration is trivial. One can then perform the standard change of variable 

{l + tf 



b = 



t 



which for < b < 4 is inverted as 



Vb-iV4~b 



Vb + i\/4-b ' 
with t varying in the unit circle between t = — 1 at 6 = and t = c~ 2l7r / 3 at b = 1; one finds 

1 e -2i^/3 

fdujf> 1] {u) = l - I *(21n(l + t)-lnt) 

-1 



-Li 2 (-t)--ln 2 £ 



Cl 2 



where use has been made of the formula 



and 



Li 2 (e^) = X -i? - X -%<t> + \f + iC[ 2 (<f,), 



<p 
Cl 2 ((f>) = - /"d» In (2 sin | J 



is the Clausen function of weight 2. 

Along the same lines, one easily obtains 



du J 2 ' (u) 



db 



du 



Vb(4 ~ b \Jv- b )° \/((i + Vbr-u)(u-(i-Vby) 



(1-v^) 2 

db , Vb-iV^b 
In- 



2 7 y/b(4-b) Vb + iV^^b 



and, as expected from Eq. f8.6fl . 



du J^\u)= ChQ 



(B.2) 

(B.3) 

(B.4) 



(B.5) 



(B.6) 



(B.7) 
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The same approach gives also at once 



du J 2 (1,9) (w) 



(B. 



The case of J} (u) is more complicated; indeed, following the same procedure as for the previous integrals 
one finds 



db 



du 



duj[ l ' 9) (u) = , , 

o ^ (4 ~ 6 V +1)2 V(«-(i + V&) 2 )(«-(i-VF) a ) 



1 



db 



8-6+y/(4-b)(16-6) 

i/6(4 - 6) n 8-6- ■ N /(4-b)(16-6) ' 



(B.9) 



where the new square root V16 — 6 has appeared. If one uses at this point the change of variable Eq.s (|B.3(l . 
as a consequence of the factor \/16 — b in the argument of the logarithm a new quadratic square root, 
vl — 14i + 1 2 , appears; that square root can be in turn removed by a suitable change of variable, and the 
result can be expressed, by brute force, as a combination of many dilogarithms of unusual and rather nasty 
arguments. It can be more convenient to consider the auxiliary integral 



A(m) 



db 



(2+m) 2 



du 



Taking the derivative with respect of to 2 of the integral representation for A(m) gives 



(B.10) 



dm 



T A(m) 



db 



2m 2 J ^6(6-4(m+l) 2 ) 



2m- 



■(—2 lnx — i7r) , 



where x is defined through m + 1 = 2x/{x 2 + 1). The primitive in x of that expression is 



A(m(x)) 



2i7rln(l - x) - «7iTn(l + x 2 ) +41na;ln(l - x) +4Li 2 (a;) 



1 



- \n{x 2 ) ln(l + x 2 ) - Li 2 (-.x 2 ) + -it 2 



(B.ll) 



where the constant term 7r 2 /4 is fixed by observing that A(m = —1) is equal to the integral Eq. (|B.7p times 
i. At to = 1, or x — cT™! 3 , Eq. (jB.ll|) gives 



Jdu4 1 *\u) = 5Cl 2 (^ 



(B.12) 



The same result can be obtained by using the identity Eq. (|9.1U|I ; by the change of variable u = 9/v one 
has 



:. 9 

du4 l ' 9 \u) = ^ dv-4 l ' 9 \v), 
9 ./ v 



and the integration can be continued, as for Ea. (|B.8Jl . along the line followed for Ea.s (|B.2IB.7p . 



(B.13) 



•SO 



With a suitable combinations of the above described techniques, one obtains also, with U » 1 

u 



duJ^°°\u) = 7rln[/-5Cl 2 (|) 

9 
U 

JduJ?'°°\u) = hr?U-\^. (B.14) 

9 

The extension of the above formulae to the case in which a logarithmic factor is also present in the 
integrand requires much additional effort. We list the results 

l 

j du J^i^lnu = 3Ls 3 (jf\ + Itt 3 , 

o 
i 

du 4 0A) (u)lnu = 2ttC1 2 (|) - 7C(3) , 



3 
o 
l 

du jf^fuMnn-w) = -— 7T 3 , 

36 

o 

i 

d«j( 0ll) («)ln(l-«) = -§C(3) 

o 

i 



JduJ^WMQ-u) = 3Ls 3 ^+^7r 3 , 



1 

J duJ^ 1] (u)H$-u) = |*C(3), (B.15) 



l 

9 



du j[^\u)\nu = 151*3 f^+H 



2tt\ 23 , 
7r d 



2 
1 
9 

(1,9) 



du J 2 (1,9) (u)lnu = 7C(3) 



du Jf '"'(«) ln(«-l) = /3 3 , 



du4 h9) (u)\n(u-l) = ^C(3), 
l 

9 

*.jf 1,9) (u)lii(9-u) = l5LsJ^]+^7r 3 +f3 3 
l 

9 



f du 4 h9 \u) ln(9 - «) = 5^C1 2 (^) - yC(3) , (B.16) 
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9 
U 



[duj[ 9 '°°\u)lnu = i7rln 2 [/-y7r 3 -15Ls 3 (y N ) , 

9 
U 

du4 9 '°°\u)lnu = iln 3 [/-C(3), 
JduJ^°°\u)ln(u-l) = K In 2 U + l^ 3 - /? 3 , 
'du4 9 -°°\u)Hu-l) = Iln 3 C/-^C(3), 

r 

Jduj[ 9 >°°\u)Mu-9) = ^ln*U-^n 3 -l5Ls 3 (jf^-(3 3 , 

r 

du4 9 -°°\u)Hu-9) = Iln 3 C/ + ^C(3)-57rCl 2 (^) . (B.17) 

The constants LS3 (^-) and j3 3 appearing in the above equations are defined as 


Ls 3 (» = - d9\n (2sin' 



9 
u 



9 
u 



2/ ' 
03 = - J^ln 2 3 - ^Li 2 f-i ) - 18 7mSi 2 (zV3) . i BJ «) 

As a final remark, in obtaining the previous results we obtained also the relations 

Ls 3 (</>)+Ls 3 (tt-</>) = -\ir\ (B.19) 

6 

where G is the Catalan constant, which we could not find in the literature. 
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